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Abstract 



Using chiral Ward identities, we determine the renormalisation constants 
of bilinear quark operators for the Sheikholeslami-Wohlert action lattice at 
/3 = 6.2. The results are obtained with a high degree of accuracy. For the vec- 
tor current renormalisation constant we obtain Zy = 0.817 it 2 it 8, where the 
first error is statistical and the second is due to mass dependence of Zy. This is 
close to the perturbative value of 0.83. For the axial current renormalisation 
constant we obtain Z^ = 1.045 + ^^, significantly higher than the value ob- 
tained in perturbation theory. This is shown to reduce the difference between 
lattice estimates and the experimental values for the pseudoscalar meson decay 
constants, but a significant discrepancy remains. The ratio of pseudoscalar to 
scalar renormalisation constants, Zp/Zg, is less well determined, but seems to 
be slightly lower than the perturbative value. 
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I. INTRODUCTION 



In a recent paper [1], we reported values for the decay constants and fx from lattice 
calculations that were considerably lower than the experimental values. This has been a 
persistent feature of lattice calculations of pseudoscalar decay constants [2,3]. However, it 
has not been possible to attribute this discrepancy to any particular systematic error in the 
calculations (eg quenching), because of uncertainties in the renormalisation of the current 
operators involved. 

A fully non-perturbative determination of renormalisation constants for finite operators 
can be achieved using chiral Ward identities [4,5], thereby bypassing the need for perturbative 
calculation of these quantities. A preliminary calculation at f3 = 6.0, presented in [5], 
indicated that the axial vector current renormalisation, in particular, differs significantly 
from its perturbative value. In this paper, we present results for the vector, axial- vector, 
pseudoscalar and scalar renormalisation constants at /5 = 6.2, and update the values for 
decay constants given in [1]. 



II. LATTICE WARD IDENTITIES AND RENORMALISATION CONSTANTS 

We define the (non-conserved) lattice vector current, the axial current and the pseu- 
doscalar and scalar densities as follows: 

(1) 
(2) 
(3) 
(4) 

In the continuum limit, these operators are related to operators obeying the correct 
current algebra by multiplicative renormalisation constants Zy, Z\ etc, so that = ZyV^ 
etc. [5]. 

The renormalisation constants for can easily be determined by evaluating 
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Inserting a complete set of states and noting that the matrix element of V^^(O) between a 
degenerate pseudoscalar meson state is (P, I F/- I P,) = 2Er,/Z{; , we see that this should 
give a precise estimate provided the effect of the off-diagonal matrix elements (P^ | \ Pn) 
can be neglected. 

For the axial case, there is no conserved current, or any other "easy" way of determining 
the renormalisation constants, but they can be obtained using chiral Ward Identities. Using 
the arguments of [4,5], we obtain the following identities for the Sheikholeslami-Wohlert 
(SW) [6] action: 
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and d"''' is defined by 
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A",A^| = 2(^"^^A^ + -^^"^ (9) 



III. COMPUTATION AND RESULTS 

We have performed simulations at /5 = 6.2 on a 24'^ X 48 lattice. We have generated 
propagators using the SW action, with two quark masses, corresponding to /c = 0.14144 
and 0.14262, where n = l/2(mo + 4r) and r = 1. 60 configurations have been analysed 
at /c = 0.14144, and 26 configurations at /c = 0.14262 (for details, see [7]). The statistical 
errors are calculated with a bootstrap procedure, using 100 bootstrap samples. 

Zv was determined from eq. (5), using 10 configurations, at three values for the quark 
mass (corresponding to /c = 0.14144, n = 0.14226 and n = 0.14262). The results are 
presented as a function of t in fig.l. We see that the values for Zy are roughly independent 
of t. Our best values, obtained by fitting to timeslices 5-19, are given in table I. The errors 
from the variation between the timeslices are obtained from fits to 100 bootstrap samples 
of timeslices within the fit range. 

The results are plotted as a function of the square of the mass of the pseudoscalar 
meson (proportional to the quark mass) in fig. 2. We see that the results show a clear 
(linear) dependence on the quark mass, consistent with the expectation that the leading 
corrections to our calculations should be of O^asirLoo). Perturbation theory at one-loop 
level [8,9] with a "boosted" coupling constant [10] g 0.83, which is quite close to 

our non-perturbative values. 
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FIG. 1. as a function of t. 
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FIG. 2. Zy as a function of the mass of the pseudoscalar meson. The scale is taken from the 
string tension [7]. 
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™P5 




0.14262 
0.14226 
0.14144 


0.208 
0.341 
0.663 


0.82139+41 + 25 
0.82453 + 24 + 24 

0.83136^23 + 22 



TABLE I. Values of the renormalisation constant Zy as a function of the quark mass. The 
first set of errors are the statistical errors, while the second set are the errors due to the variation 
between the timeslices. 
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FIG. 3. as a function of i for k = 0.14144 
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Zp/Zs 


0.14262 
0.14144 


0.208 
0.663 


1.040 + 1° 
1.047+^ 


0.693128 
0.6491^ 



TABLE II. Values of the renormalisation constants Z^ and ZpjZg as functions of the quark 
mass. 



We have also used eq. (6), with = /) = 0, to check the consistency of our results at 
K = 0.14144, using the value for Zy quoted in table I and the value for given in table II 
as input. This gives Zy = 0.817 1^^, which is within 2a of the result obtained from eq. (5). 

The axial vector renormalisation constant Z^ is determined using eq. (6), with u = p = i 
and summing over i = 1, 2, 3, using the values for Zy quoted in table I as input. The results 
for K, = 0.14144 are plotted against t in fig. 3. We see that, apart from the effect of the 
contact terms on the first few timeslices, they show virtually no dependence on t. Our best 
estimates are given in table II. 

Within the statistical errors, these results show no dependence of Za on the quark 
mass. The comparison with results from perturbation theory is more interesting: one-loop 
calculations with a "boosted" coupling constant g 0.97, which is considerably lower 

than our non-perturbative results. The discrepancy is higher at lower /5, as expected; in [5] 
Za at /5 = 6.0 was found to be 1.09. 

In table III we show how the values for the decay constants reported in [1] change 
when we use the results given above for the renormalisation constants. For Zy, we have 
extrapolated the values in table I to the limit of zero quark mass, giving Zy = 0.817 it 2, 
with an additional uncertainty due to the quark mass dependence of Zy of ±0.008, which 
corresponds the difference between the value at our largest quark mass and the value for 
zero quark mass. For Za, we have taken a best estimate, combining our results at the two 
/c- values, of Za = 1.045 ^1°, with the errors corresponding to the spread between the highest 
and lowest estimate. We see that all the decay constants move closer to the experimental 
values, but that a significant discrepancy still remains, especially for and fx- f-w turns 
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0.172 




0.160 + ^ 


0.172 + 8 


0.208 




0.144+4 


0.155+5 


0.179 



TABLE III. Values of decay constants in physical units, using perturbative and non-perturbative 
values for the renormalisation constants. The second set of errors in the vector meson decay con- 
stants are systematic uncertainties due to the quark mass dependence of Zy. 



out to be about 3(7 away from its experimental value. The APE collaboration has found 
f-w I {rripZA) = 0.186(20) at /5 = 6.2 [11], which gives a value for f-w/rrip compatible with 
experiment. 

The ratio of pseudoscalar to scalar renormalisation constant is determined using (7). 
The results are given in table II. The uncertainty in these results is too large to determine 
whether there is any dependence on the quark mass here. Perturbative calculations with a 
"boosted" coupling constant give Zp/Zs = 0.68. As can be seen, our result for the heavier 
quark mass (which is the more accurate) is slightly lower than this, while the lighter quark 
mass gives a value compatible with perturbative results (although the errors here are still 
quite large). Comparison with the result reported in [5] at /5 = 6.0 shows that, as in the 
case with Za, the discrepancy decreases with increasing f3. 

IV. CONCLUSIONS 

In this paper we have reported the determination of lattice renormalisation constants us- 
ing chiral Ward identities with the Sheikholeslami-Wohlert action. Our results are obtained 
with good accuracy, yielding values for Zy close to the values from perturbation theory 
(but increasing with increasing quark mass), while the value for Za is considerably higher 
than perturbative results. For Zp/Zs the uncertainties are larger, but the results we can 
have confidence in lie slightly lower than perturbative values. The result for Za brings our 
estimates for and fx considerably closer to experimental values — within 3(7 for 
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